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We study the mass of the state Y~(2175) of J PC = 1 — _ in the QCD sum rule. We construct 
both the diquark-antidiquark currents (ss)(ss) and the meson- meson currents (ss)(ss). We find 
that there are two independent currents for both cases, and derive the relations between them. The 
OPE convergence of these two currents is sufficiently fast, which enables us to perform good sum 
rule analysis. Both the SVZ sum rule and the finite energy sum rule lead to a mass around 2.3 ±0.4 
GeV, which is consistent with the observed mass within the uncertainties of the present QCD sum 
rule. The coupling of the four-quark currents to lower lying states such as 0(1020) turns out to be 
rather small. We also discuss possible decay properties of Y(2175) if it is a tetraquark state. 
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O I. INTRODUCTION 

m ' 

The theory of the strong interactions, Quantum Chromodynamics (QCD), originated from the systematics of hadron 
spectroscopy. The spectroscopy contains meson and baryon states, many of which are well classified by the quark 
Cln' model with the quark content qq and qqq. Besides the quark model, QCD allows much richer hadron spectrum such 
yi^ , as multiquark states, hadron molecules, hybrid states, glueballs etc. However the spectrum of QCD seems to saturate 
' at qq and qqq. Since 2003, there have been important developments in hadron spectroscopy, which is triggered by the 
observation of the pentaquark + . After three years of intense study, the status of + is still controversial. However, 

1 we have the charm-strange mesons £> sJ (2317), D sJ (2460) (JIU; the charmonium state Y(3872) Q. Y(4260) Q, and 

I many Xs and Ys, whose properties seem difficult to be explained by the conventional picture of qq. 

Recently Babar Collaboration observed a resonance Y(2175) near the threshold in the process e + e~ — > 0/o(98O) via 
£C) ' initial-state radiation It has the quantum numbers J PC = \~~ . The Breit-Wigner mass is M = 2.175±0.010± 

0.015 GeV, and width is V — 0.058±0.016±0.020 GeV. It has been also confirmed by BES collaboration in the process 
J/yb -> #/ (980). A fit with a Breit-Wigner function gives the peak mass and width of M = 2.186 ± 0.010 ± 0.006 
^ i GeV and T = 0.065 ± 0.023 ± 0.017 GeV Q. 

t-H ■ There are many suggestions to interpret this resonance. Ding and Yan interpreted it as a strangeonium hybrid 
and studied its decay properties in the flux-tube model and the constituent gluon model. Furthermore, for testing 
ssg scenario, they suggested searching decay modes such as Y(2175) — * Ki(lAQQ)K — > nK*(892)K, Y(2175) — > 
K 1 {1270)K -> pKK and Y(2175) -► K 1 (1270)K -> 7^(1430) @. In Ref. [U%, the authors explored T(2175) as a 
2 3 £>i ss meson, and calculated its decay modes by using both the 3 Po model and the flux-tube model. They suggested 
experimental search of the decay modes KK, K*K*, K(1460)K and hi(l380)i]. The characteristic decay modes of 
V(2175) as either a hybrid state or an ss state are quite different, which may be used to distinguish the hybrid and ss 
schemes. Wang studied T(2175) as a tetraquark state ssss by using QCD sum rule and suggested that there may be 
some tetraquark components in the state V(2175) [ll|. In a recent article [13], Zhu reviewed Y(2175) and indicated 
that the possibility of Y(2175) arising from S-wave threshold effects can not be excluded. Napsuciale, Oset, Sasaki and 
Vaquera-Araujo studied the reaction e + e~ — > (jnnr for pions in an isoscalar S-wave channel which is dominated by the 
loop mechanism. By selecting the </>/o(980) contribution as a function of the e + e~ energy, they also reproduced the 
experimental data except for the narrow peak [131 ]. Bystritskiy, Volkov, Kuraev, Bartos and Secansky calculated the 
total probability and the differential cross section of the process e + e~ — > 0/o(98O) by using the local NJL model 16 1. 
Anikin, Pire and Teryaev studied the reaction 7*7 — > pp, and calculated the mass of the isotensor exotic meson 
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In Ref. [151 ] . the authors performed a QCD sum rule study for 1 hybrid meson, and the mass is predicted to be 
2.3 — 2.4, 2.3 — 2.5, and 2.5 — 2.6 GeV for qqg, qsg, and ssg, respectively. 
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In this work, we revisit the possibility of Y(2175) as an tetraquark state ssss. With the approach developed in our 
previous work [171 ] ■ we construct the general tetraquark interpolating currents with the quantum numbers J = 1 . 
We find that there are two independent currents. They can have a structure of diquark-antidiquark (ss)(ss), or have 
a structure of meson-meson (ss)(ss). We show that they are equivalent, and derive the relations between them. Then 
by using these two independent currents, we also perform a QCD sum rule analysis. We calculate the OPE up to the 
dimension 12, which contains the (qq) condensates. In these two respects, our study differs from the previous one of 
Ref. [IH. 

This paper is organized as follows. In Sec. [Ill we construct the tetraquark currents using both diquark (qq) and 
antidiquark (qq) fields, as well as quark-antiquark (qq) pairs. In Sec. IIII[ we perform a QCD sum rule analysis by 
using these currents. In Sec. IIV1 the numerical result is obtained for the mass of Y(2175). In Sec. \V\ we use finite 
energy sum rule to calculate its mass again. Sec. IVIl is a summary. 



II. INTERPOLATING CURRENTS 



In this section, we construct currents for the state Y(2175) of J PC = From the decay pattern Y(2175) — > 

</>(1020)/o(980), we expect that there is a large ssss component in Y(2175). We may add further quark and antiquark 
pairs, but the simplest choice would be ssss. We will discuss later how this simplest quark content will be compatible 
with the above decay pattern when considering the possible structure of 0(1020) and /q(980). 

Let us now briefly see the flavor structure of the current. In the diquark-antidiquark construction (ss)(ss) where 
ss and ss pairs have a symmetric flavor structure, the flavor decomposition goes as 

6 f ® 6 f = If 9 8 f © 27 f . (1) 

Therefore, the (ss)(ss) state is a mixing of l/,8/ and 27/ multiplets in the ideal mixing scheme. 

Now we find that there are two non- vanishing currents for each state with the quantum number J PC = 1 . For 
the state ssss: 

Vin = { s lC^ 5 s b )(s allll5 Csl) - (slC^^r oSb )(s al5 Csl) , (2) 
V2n = {slC^Sh^Saa^Cs^) - (s'lC(J llv s b )(s a 'y u Cs b r ) , (3) 

where the sum over repeated indices (/i for Dirac spinor indices, and a, b for color indices) is taken. C — 17270 is the 
Dirac field charge conjugation operator, and the superscript T represents the transpose of the Dirac indices only. 

Besides the diquark-antidiquark currents, we can also construct the tetraquark currents by using quark-antiquark 
(ss) pairs. We find that there are four non-vanishing currents: 

»73m = (SaSa)(sbl^Sb) , 

Vip = (Sa7 l/ 75Sa)(s&Cr Al „75S6) , 

Vs? = A ab A cci (s a s & )(s c 7 A1 s (i ) , 

V6fj, = KbKd(Sal U lbSb)(s c CrpvlbSd) ■ 

In Ref. pj| . the author used 775^, to perform QCD sum rule analysis, which is a mixing of 771^ and 772^. We can verify 
the following relations by using the Fierz transformation: 

»/6p = - , mn = 3z7 73 M + ■ (4) 

Therefore, among the four (qq)(qq) currents, two are independent. We can also verify the relations between (ss)(ss) 
currents and (ss)(ss) currents, by using the Fierz transformation: 

"Hin = -V3n + , rj2» = 3«773 M - V4^ ■ (5) 
Therefore, these two constructions are equivalent, and we will use 771^ and 772^ for QCD sum rule analysis. 



III. QCD SUM RULE ANALYSIS 



For the past decades QCD sum rule has proven to be a very powerful and successful non-perturbative method [H,[l||, 
and it has been applied to study tetraquark states in many references [13; US HH, HI, HI, [24, 25]. In sum rule analyses, 
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we consider two-point correlation functions: 

n^( q 2 ) = i J ^^(01^(^(0)10) , (6) 

where rj^ is an interpolating current for the tetraquark. The Lorentz structure can be simplified to be: 

IW? 2 ) = - ^)n«(<z 2 ) + M^n(°V) . (7) 
i i 

We compute H(q 2 ) in the operator product expansion (OPE) of QCD up to certain order in the expansion, which 
is then matched with a hadronic parametrization to extract information of hadron properties. At the hadron level, 
we express the correlation function in the form of the dispersion relation with a spectral function: 

n (1) (g 2 )=f P( 2 S) . da, (8) 

where the subscript is (4m s ) 2 = 16m 2 , and 

p(s) ee ^5( a -M*)(0\ri\n){n\rf\0) 

■n 

= fy8(s - My) + higher states . (9) 

For the second equation, as usual, we adopt a parametrization of one pole dominance for the ground state Y and a 
continuum contribution. The sum rule analysis is then performed after the Borel transformation of the two expressions 
of the correlation function, © and © 

/•OO 

n(«») ( Ml) = B MS n( 1 )(p 2 )= / e-°/ M Bp{s)d S . (10) 

Assuming the contribution from the continuum states can be approximated well by the spectral density of OPE above 
a threshold value sq (duality), we arrive at the sum rule equation 

n(M|) = /2 e -A4/M£ = / e s/Ml p{s)ds _ (n) 

Jl6m 2 s 

Differentiating Eq. (fTTj) with respect to 1/M B and dividing it by Eq. (fTTj) . finally we obtain 

n(Afl) j^e-^lp{s)ds ■ [12) 

In the following, we study both Eqs. (fTTj) and (fT2|) as functions of the parameters such as the Borel mass Mb and the 
threshold value sq for various combinations of the tetraquark currents. 

For the currents rji^ and 772^, we have calculated the OPE up to dimension twelve, which contains the {qq) A 
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condensate: 

ni (Mi) 
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(13) 
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(14) 



In the above equations, (ss) is the dimension D = 3 strange quark condensate; (g 2 GG) is a D = 4 gluon condensate; 
(gsaGs) is -D = 5 mixed condensate. There are many terms which give minor contributions, such as (g 3 G 3 ), and 
we omit them. As usual, we assume the vacuum saturation for higher dimensional condensates such as (0\qqqq\0) ~ 
(0|g<7|0)(0|gg|0). To obtain these results, we keep the terms of order 0(m 2 ) in the propagators of a massive quark in 



the presence of quark and gluon condensates: 

iS ah ee (0\T[q a (x)q\0)]\0) 



id 



a b 



A 



ah 



5 ab 



27T 2 5 



5 ab x 2 



32tt 2 2 
(g c qaGq) - 



m q 5 ah iS ab m q (qq) 



12 

ib ah m\ 



192 " 47r 2 a: 2 48 8tt 2 x 2 

We find that there is an approximate relation between the correlation functions of rji^ and 772^: 

3ni(Af|)~2n 2 (M|), 



(15) 



(16) 



which is valid for the continuum, (ss), and (g c qaGq) terms, etc. So the numerical results by using them are also very 
similar. 



IV. NUMERICAL ANALYSIS 



In our numerical analysis, we use the following values for various condensates and m s at 1 GeV and a s at 1.7 GeV 

m si m, m, m, m m hi : 

(qq) = -(0.240 GeV) 3 , 

(ss) = -(0.8 ± 0.1) x (0.240 GeV) 3 , 

(g 2 GG) = (0.48 ± 0.14) GeV 4 , 

( 9s qaGq) = -M 2 x (qq) , (17) 
M 2 = (0.8 ± 0.2) GeV 2 , 
m s (l GeV) = 125 ± 20 MeV, 
a s (1.7GeV) = 0.328 ± 0.03 ± 0.025 . 
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There is a minus sign in the definition of the mixed condensate (g s qaGq), which is different from that used in some 
other QCD sum rule studies. This difference just comes from the definition of coupling constant g s [27l. 1341]. 

First we want to study the convergence of the operator product expansion, which is the cornerstone of the reliable 
QCD sum rule analysis. By taking sq to be oo and the integral subscript f 0>m 2 s to be zero, we obtain the numerical 
series of the OPE as a function of M B '■ 

n x (Mf ) = 1.4 x IQ- & M B ° - 3.8 x 10~ 7 A/f - 6.2 x 10~ 7 Aff + 4.2 x 10" 7 M B 

-1.2 x lCT 6 AfJ> + 4.7 x 1(T 8 - 1.5 x IQ- 7 M B 2 , (18) 

n 2 (Af|,) = 2.0 x 1(T 6 M^ - 5.7 x 1CT 7 ' M% - 8.0 x 1Q- J M B + 6.4 x 1Q- 7 M| 

-1.7 x 10~ 6 A/f + 1.0 x 10~ 7 - 2.2 x 10~ 7 M B 2 . (19) 

After careful testing of the free parameter Borel mass M B , we find for AfJ, > 2 GeV 2 , which is the region suitable 
for the study of V(2175), the Borel mass dependence is weak. Moreover, the convergence of the OPE is satisfied in 
this region. The correlation function of the current ryi„ is shown in Fig. [TJ when we take Sq — 5.7 GeV 2 (the integral 
subscript is still 16m 2 ). We find that in the region of 2 GeV 2 < M B < 5 GeV 2 , the perturbative term (the solid line 
in Fig. [T|) gives the most important contribution, and the convergence is quite good. 




FIG. 1: Various contribution to the correlation function for the current 77i M as functions of the Borel mass Mb in units of 
GeV 10 at so = 5.7 GeV 2 . The labels indicate the dimension up to which the OPE terms are included. 



It is important to note that the Y(2175) state is not the lowest state in the 1~~ channel containing ss and that 
the interpolating currents see only the quantum number of the states. It is possible that the lowdying states 0(1020) 
and 0(1680) also couple to the tetraquark currents 771^ and 772^. If so, their contribution to the spectral density and 
the resulting correlation function should be positive definite. 

However, we find that (1) the spectral densities p(s) of Eq. §§§ for both currents rji^ and 7/2^ are negative when 
s < 2 GeV 2 ; (2) the Borel transformed correlation function U(M B ) in Eq. (fTTj) is also negative in the region so < 4.3 
GeV 2 and 1 GeV 2 < M B < 4 GeV 2 . As an illustration, we show the correlation function as a function of so in Fig. [5] 
This fact indicates that the ssss tetraquark currents couple weakly to the lower states 0(1020) and 0(1680) in the 
present QCD sum rule analysis. 

The pole contribution is not large enough for both currents due to D = 10 perturbative term e~ s l M B s 4 ds, which 
is a common feature for any multiquark interpolating currents with high dimensions. The mixing of the currents 7/1^ 
and r]2fj, does not improve the rate of the pole contribution. The small pole contribution suggests that the continuum 
contribution to the spectral density is dominant, which demands a very careful choice of the parameters of the QCD 
sum rule. In our numerical analysis, we require the extracted mass have a dual minimum dependence on both the 
Borel parameter M B and the threshold parameter sq. In this way, we can find a good working region of M B and sq 
(Borel window), where the mass of V(2175) can be determined reliably. 

Now the mass is shown as functions of the Borel mass M B and the threshold value sq in Fig. [3] and Fig. |4j The 
threshold value is taken to be around 5^7 GeV 2 , where its square root is around 2.2 ~ 2.7 GeV. We find that there 
is a mass minimum around 2.4 GeV for the current 771^, when we take M B ~ 4 GeV 2 and so ~ 5.7 GeV 2 . While this 
minimum is around 2.3 GeV for the current 772^, when we take M B ~ 4 GeV 2 and so ~ 5.4 GeV 2 . 

In short summary, we have performed the QCD sum rule analysis for both 771^ and 772^. The obtained results 
are quite similar. This is due to the similarity of the two correlation functions as shown in Eq. (|16[) . We have also 
considered their mixing, which also give the similar result. The mass is predicted to be around 2.3 ~ 2.4 GeV in the 
QCD sum rule. 
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FIG. 2: The correlation function for the current 771^, as a function of so in units of GeV 10 . The curves are obtained by setting 
M% = 1 GeV 2 (long-dashed line), 2 GeV 2 (short-dashed line) and 4 GeV 2 (solid line). 




s [GeV^ ] Borel Mass [GeV^ ] 

FIG. 3: The mass of Y~(2175) as a function of Mb (Left) and so (Right) in units of GeV for the current 

V. FINITE ENERGY SUM RULE 

In this section, we use the method of finite energy sum rule (FESR). In order to calculate the mass in the FESR, 
we first define the nth moment by using the spectral function p(s) in Eq. ([9]) 



W(n, s ) 



p(s)s n ds . 



(20) 



This integral is used for the phenomenological side, while the integral along the circular contour of radius Sq on the 
q 2 complex plain should be performed for the theoretical side. The lower integral bound s = is taken in order to 
include the deZta-function contribution in the OPE (Eqs. (|2"5)l and l|24p). 
With the assumption of quark-hadron duality, we obtain 



W(n, s ) 

The mass of the ground state can be obtained as 



Hadron 



W(n, s ) 



OPE 



M^(n,s ) 



W(n+l,s ) 
W(n,s ) 



(21) 



(22) 



For the currents 771^ and 772^, the spectral functions pi(s) and P2(s) can be drawn from Eqs. (fT5)) and (fTij) . The d 
= 12 terms which are proportional to l/(q 2 ) 2 do not contribute to the function W(n,so) of Eq. (f2"0"|) for n = 0, or 
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FIG. 4: The mass of Y(2175) as a function of Mb (Left) and so (Right) in units of GeV for r\2^- 



they have a very small contribution for n = 1, when the theoretical side is computed by the integral over the circle 
of radius sq on the complex q 2 plain. Therefore, the spectral densities for 771^ and 772^ take the following form up to 
dimension 10, 
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FIG. 5: The mass of F(2175) by using the current 77i M (solid line) and 772^ (dashed line) as a function of so in units of GeV. 



The mass is shown as a function of the threshold value so in Fig. [SJ where n is chosen to be 1. We find that there 
is a mass minimum. It is around 2.3 GeV for the current rji^ when we take sq ~ 5.2 GeV 2 , while it is around 2.2 
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GeV for the current r/2 M when we take sq ~ 4.8 GeV 2 . For the current 771^, the minimum point occurs at t/sq = 2.28 
GeV where the mass takes 2.3 GeV, and the threshold value is slightly smaller than the mass, unlike the ordinary 
expectation that t/sq is larger than the obtained mass. However, the minimum point is on the very shallow minimum 
curve and the resulting mass is rather insensitive to the change in the y/Jo value. Therefore, we can increase y/s~o 
slightly more, for example 2.45 GeV, but the mass still remains at around 2.35 GeV, which is smaller than ^/sq now. 
This fact is due to the uncertainty of our sum rule analysis as well as the negative part of the spectral densities. For 
example, if we take the lower limit of integrations in Eq. (|2"2l to be 1 GeV 2 instead of 16m 2 (the positive part starts 
at around 3 GeV 2 ), the mass minimum will be around 2.1 GeV, when Sq is around 4.5 GeV 2 ; if we take the lower limit 
to be 2 GeV 2 , the mass minimum would be around 2.0 GeV, when sq is around 4 GeV 2 . We show the second case in 
Fig. M The region 5 < s < 6 GeV 2 is suitable for the QCD sum rule analysis, and the mass obtained is around 2.2 
GeV. Therefore, considering the uncertainty of the QCD sum rule, we obtain the same result as the previous one. 

3.0 



% 2.5 
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% 2.0 
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3 4 5 6 
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FIG. 6: The mass of V(2175) by using the current 771^ (solid line) and 772^ (dashed line) as a function of so in units of GeV, 
when the lower limit of integrations in Eq. (122[l is 2 GeV 2 instead of 16m 2 . 



VI. SUMMARY AND DISCUSSIONS 

In this work we have studied the mass of the state V(2175) with the quantum numbers J PC = 1~~ in the QCD sum 
rule. We have constructed both the diquark-antidiquark currents (ss)(ss) and the meson-meson currents (ss)(ss). 
We find that there are two independent currents for both cases and verify the relations between them. Then using 
the two (ss)(ss) currents, we calculate the OPE up to dimension twelve, which contains the (ss) 4 condensates. The 
convergence of the OPE turns our to be very good. We find that the OPE's of the two currents are similar, and 
therefore, the obtained results are also similar. By using both the SVZ sum rule and the finite energy sum rule, we 
find that there is a mass minimum. For SVZ sum rule, the minimum is in the region 5 < sq < 7 GeV 2 and 2 < AfJ, < 4 
GeV 2 . For finite energy sum rule, the minimum is in the region 4.5 < sq < 5.5 GeV 2 . Considering the uncertainty, 
the mass obtained is around 2.3 ±0.4 GeV. The state Y(2175) can be accommodated in the QCD sum rule formalism 
although the central value of the mass is about 100 MeV higher than the experimental value. We calculate the OPE 
up to dimension twelve and include many terms, but still the accuracy is around 20%. This is the usual accuracy of 
the QCD sum rule. In our analysis it is partly due to the many omitted condensates such as (GGG) etc. 

We have investigated the coupling of the currents to the lower lying states including 0(1020) and found that the 
relevant spectral density becomes negative, implying that the present four-quark currents can not describe those states 
properly. This fact indicates that the four-quark interpolating currents couple rather weakly to 0(1020), which is a 
pure ss state. 

We can test the tetraquark structure of Y(2175) by considering its decay properties. Naively, the ssss tetraquark 
would fall apart via S'-wave into the 0(lO2O)/o(98O) pair, and would have a very large width. The experimental width 
of Y(2175) is only about 60 MeV, which seems too narrow to be a pure tetraquark state. We can discuss the decay 
of the y(2175) by borrowing an argument based on a valence quark picture. The (ss)(ss) configuration for V(2175) 
can be a combination of 3 Si and 3 Poi which may fall apart into two mesons of 1~ and + in the s-wave. In the QCD 
sum rule the 1~ ss meson is well identified with <^(1020) , while the + ss meson has a mass around 1.5 GeV and is 
hard to be identified with the observed /o(980). Therefore, such a fall-apart decay would simply be suppressed due to 
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the kinematical reason. The physical /o(980) state may be a tetraquark state as discussed in the previous QCD sum 
rule study [17l] . Then the transition Y(2175) — > 0(1020) + f (tetraquark) should be accompanied by a qq creation 
violating the OZI rule, as well as by an annihilation of one quanta of orbital angular momentum. These facts may 
once again suppress the decay of Y(2175) — > 0(1020) + /o(980). This fact was studied in the recent paper by Torres, 
Khcmchandani, Geng, Napsuciale and Oset [35|. They studied the tf>KK system with the Faddeev equations where 
the contained KK form the /o(980) resonance. The decay width they calculated is around 18 MeV, not far from the 
experimental value. The all above evidences would imply that the Y(2175) is a possible candidate of a tetraquark 
state. 

Y(2175) could be a threshold effect, a hybrid state ssG, a tetraquark, an excited ss state or a mixture of all the 
above possibilities. Because of its non-exotic quantum number, it is not easy to establish its underlying structure. 
Clearly more experimental and theoretical investigations are required. 

One byproduct of the present work is the interesting observation that some type of four-quark interpolating currents 
may couple weakly to the conventional qq ground states. If future work confirms this point, we may have a novel 
framework to study the excited qq mesons using the four-quark interpolating currents, which is not feasible for the 
traditional qq interpolating currents. 
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